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ABSTRACT 

We revisit the shocked shell model for the class of Active Galactic Nuclei 
known as Gigahertz Peak Spectrum sources, incorporating new observational 
data on the radiation brightness temperatures. We argue that in addition to 
free-free absorption, induced Compton scattering will also have an important 
effect in forming the ~ GHz peak and in shaping the radio spectra that 
characterize these sources. Indeed, our arguments suggest that GPS sources 
may provide the first real evidence for the role of induced Compton scattering 
in extragalactic radio sources. 

Subject headings: galaxies: active — shocks — radiation mechanisms: thermal, 
non-thermal — scattering 

1. Introduction 

Recently, Bicknell, Dopita, & O'Dea (1997) developed a model which unifies the 
Gigahertz Peak Spectrum (GPS) sources and the related Active Galactic Nuclei (AGN) 
classes known as Compact Steep Spectrum (CSS) sources and Compact Symmetric Objects 
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(CSOs). Their model attributes the observed radio spectra and optical emission to the 
interaction of a jet-driven, non-thermal lobe with the ambient galactic interstellar medium 
(ISM). As the radio lobe forces its way through the ISM, strong radiative shocks create 
shells of shock-ionized and photo-ionized gas which produce a low-frequency break 
in the power law spectrum through free-free absorption (FFA). Although synchrotron 
self-absorption (SSA) provides a natural mechanism to produce a low frequency break in 
a non-thermal spectrum and is successful in explaining the spectra of AGN with distinct 
core-jet morphologies ( piandford fc Konigl 1979| ), we expect that it is unimportant in the 
more extended jet lobes of GPS sources, where the non-thermal plasma is more likely to be 
optically-thin to its own radiation. However, the spectrum emitted by the non-thermal lobe 
plasma can still exhibit a low frequency break not only as a result of FFA in an external 
shell of thermal gas surrounding the lobe, as suggested by Bicknell et al. (1997), but also 
as a result of the competing low-frequency process of induced Compton scattering (ICS). 

ICS becomes important in a high brightness temperature radiation field, such as 
that which is readily generated by synchrotron radio sources and which has a Thomson 
scattering optical depth and a brightness temperature kT^^ i^ ^ mgC^ 3> hv, such that 

f/^)r.^l, (1) 
V rricC^ ) 

where the parameter / quantifies the strong effect of the angular distribution of the 
radiation and is given by: 



f^— (I + cos' <i)){l- cos <p)dVL (2) 

loTT J 

where (f) is the scattering angle and the integral is over the solid angle of the beam incident 
on the scattering region ( |Coppi, Blandford fc Rees 1993| ). For a conical beam of radiation 
with half angle cos~ l(/io) and sohd angle f2o, / = (1 ~ /^o)^(3/Uq + 2/io + 7)/32; for complete 
isotropy of the radiation field / = 1; for /xq near unity, / ~ |(1 — /io)^- In both cases 
/ (f2o/47r)^. ICS redistributes photon energies from frequencies near where T^^y peaks 
to lower frequencies where they are absorbed (e.g. by SSA or FFA) effectively reducing 
the peak T^.^^, observed. As first demonstrated by [Sunyaev (1970)1 , the key observational 



consequence is a break in synchrotron radio spectra. 

Equation (|1]) implies that induced scattering can produce spectral distortions when 
Tt Jc; (5 X 10^)K//Tb and hence, unlike 'ordinary' (spontaneous) Compton scattering, it 
can be important even when < 1. Indeed, for a typical synchrotron radio source, 
ICS can be more efficient than SSA in producing a spectral turnover and can limit the 
brightness temperature at GHz frequencies to Tb^icHz ~ 10^^t'9'°'^7j^fn where 7min is the 
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low-energy cutoff to the non-thermal electron distribution! (^incell fc Krolik: 1994|) . Note 
that this limit is only weakly dependent on the source parameters and (for / ~ 1) is lower 
than the theoretical inverse Compton limit, < 10^^ K, for typical AGN sources ( [Hoyle,| 
Burbidge fc Sargent 1966| ). As argued by Sincell & Krolik (1994), ICS may explain why the 



inverse Compton limit is, in practice, rarely observed. Such is the case for GPS sources, 
where the observed brightness temperatures of lobe subcomponents are typically 10^~^° K 
( ptanghellini et al. 1997|) . Since ICS by the non-thermal, relativistic electrons in the jet 



lobes cannot reduce the intrinsic peak to these observed levels, then ICS may instead 
be operating externally. Indeed, Coppi, Blandford & Rees (1993) have shown that ICS by 
relatively 'cold' (nonrelativistic) electrons outside of a primary non-thermal source can 
effectively reduce the peak brightness temperature to values such that kT^^ ~ mgC^ and this 
alone makes it very tempting to adopt such a scenario to explain the observed values in 
GPS sources. 

In this Letter, we examine the process of induced scattering to show that in order to 
account for the observed brightness temperatures in GPS sources, this process must be 
operating at some level in the shock-ionized shell that has been postulated to surround the 
non-thermal lobes. 



2. The Case for Induced Scattering in GPS Sources 



In the scenario proposed by Picknell et al. (1997)| , the non-thermal jet lobes in GPS 



sources are surrounded by a shell of thermal, nonrelativistic plasma formed from shocked 
and photoionized material in the ambient ISM. Although FFA is likely to be operating at 
some level in this external shell, there will be strong competition from ICS, especially if 
the shell is tenuous. Indeed, with the expected brightness temperatures, ~ 10^^ K, of the 
primary non-thermal radiation which impinges upon the shell, ICS is more efficient than 
FFA in producing a break at ~ GHz frequencies when the electron number density, Ug, 
satisfies 

ne«100.,^Tf^(^)cm-3, (3) 

where the electron temperature, Te = lO'^T^ K. This is precisely the range of densities 
relevant to shock-ionized ISM material surrounding jet lobes in GPS sources, as found 
by picknell et al. (1997)| . Although, in general, collective plasma effects will not have an 



•^A characteristic power-law index of 2.4 has been assumed here for the electron distribution. According 
to standard synchrotron theory, this corresponds to a spectral index of a ~ 0.7, which is a fiducial value for 
the GPS sources we are concerned with here (Bicknell et al. 1997). 
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overall affect on ICS (see Uoppi, Blandford fc Rees 1993| ), this is not strictly the case in 



the presence of a high brightness temperature radiation field, which can trigger stimulated 
Raman scattering, whereby the radio waves scatter off plasmons (Langmuir waves) rather 
than off the individual electrons themselves. Since, at the densities relevant here, this 
process requires brightness temperatures > 10^^ K ( [Levinson fc Blandford 1995|) , it is 
unlikely to be important in this context. 

A definitive test for any model which attempts to explain the spectral peak in GPS/CSS 
sources is the observed inverse correlation between the frequency at which the spectrum 
peaks, z/p, and the source size. In the FFA model of [Bicknell et al. (1997)| , a theoretical 
relationship between these two parameters readily emerges from the shock dynamics. We 
follow the same method here in constructing an analogous ICS model, but as we explain 
below, the predicted relationship has different dependences on the physical parameters. 

A spectral break due to ICS is expected at the frequency where the optical depth, 
given by the left hand side of equation reaches unity. However, as equations (0) and (0) 
show, induced scattering depends strongly on the angular distribution of the high brightness 
temperature radiation field and we have allowed for this through the parameter /. The flux 
density corresponding to the condition on the brightness temperature is obtained by simply 
integrating the surface brightness over the solid angle Q subtended by the source at the 
observer, giving 

i^~^F^ = 2me / (/ T^y^ dn . (4) 



In the mo del developed by |Bicknell et al. (1997)| , Tt is approximately uniform over the sides 
of the lobe and the solid angle can be expressed as Afi ~ nrcXh/D"^ = tt^^^^'^x'^/ D"^, where 
Tc is the radius of the semi-ellipsoidal cocoon which forms as the jet-driven lobe expands 
into the ISM, is the distance from the core at the base of the jet to the lobe hotspot (i.e. 
the semi-major axis of the cocoon), D is the distance to the source and C(— 2 typically) 
is the ratio of the averaged hotspot pressure to cocoon pressure. By taking a power-law 
spectrum, Fi, = Fy^^{v / v^-)^" , where is a fiducial frequency (5 GHz in this case), and using 
the monochromatic power P^^ = AttD^F^^, equation (H) predicts a peak frequency given by 



UO ' T 



1 

a + 2 



(5) 



where P^^ = AnD F^^ is the monochromatic power. 



The Thomson optical depth is assumed to be approximately uniform over the sides 
of the lobe and is calculated with the MAPPINGSII code ( Sutherland fc Dopita 1993| ) for 
radiative shocks, giving 

Tt ~ (4.9 X 10"^) Kf-^ + (1.6 X 10-2) V^-^ , (6) 
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where V = lOOOVakms"^ is the shock velocity corresponding to the sideways expansion of 
the radio lobe. The two terms correspond to the shocked gas and photoionized precursor. 
Specifically, the shock velocity is given by (see equations 2-10 and 2-11 in Bicknell et al. 
1997) 



<5-2 
3 



where 

V,, = l.6(^)"\^'^\^f (8) 

V 8 — / I no . 

and where Fe = 10^^Fe,45 ergs~^ is the jet energy fiux and 6 is the power-law index of the 
Hydrogen density distribution, given by nn = no(xh/kpc)^'^. 

Equation (|^) predicts a relationship between peak frequency and source size that 
derives quite differently from that predicted by the FFA model (see equation 5-10 in 
Bicknell et al. (1997)|) . Since the opacity of ICS is proportional to the radiation brightness 



temperature, the chief dependence of Up derives from the integrated surface brightness 
(oc x^). There is also a weaker dependence on Xh through the shock velocity, as given in 
equation (|^, which in turn determines Tt. In the case of the FFA model, on the other hand, 
while there is also a similar weak dependence through V, the predicted Up - relationship 
derives chiefly from the dependence on density (which has a power-law scaling with source 
size). This density dependence arises simply because the opacity due to FFA is cx / nidi, 
whereas for ICS, it is oc / ngdl where both integrals are through the ionized screen. Allowing 
for the inverse dependence of the shocked cooling zone and the photoionized precursor on 
density, the dependences on density are respectively on the flrst and zeroth powers of ng. 
This is reflected, for instance, in the dependence of the Thomson optical depth on velocity 
alone and this is where a dependence on density does enter, since the expansion velocity is 
a function of the ambient density. Although these two different models predict a z/p X]^ 
relationship that derives from intrinsically different physical processes (i.e. two-body 
interactions versus particle-photon interactions), the overall differences are surprisingly 
small: both FFA and ICS models predict a negative slope in the i^p - x^ plane which is 
close to unity and which varies very weakly with the relevant range of parameters. 



Let us now examine the predicted peak frequency - source size relationship more 
quantitatively. 



- 6 - 



3. Results 

The relationship between peak frequency, z/p, and source size, Xh, as predicted by 
equation for the anisotropy parameter / = 1 is plotted and compared with observational 
data ( Fanti et al. 19901 , gtanghellini et al. 199"^ , P'Dea fc Baum 19971 ) in figure 1. The 



parameters involved in the theoretical lines are as follows: values of the Hydrogen densities 
at 1 kpc, no, are 1 and lOcm"'^; values of the density power-law index, 6, are 1.5 and 2 (these 
are the lower and upper limits suggested by Begelman (1996) in his fit to the luminosity-size 



relation for CSOs); the total 5 GHz radio power is Psghz = 10^^'^ erg s ^Hz ^, corresponding 
to the average radio power of GPS sources ( [t'anti et al. 1990| , |Stanghellini et al. 1997| ); the 



jet energy flux is = 10''^ ergs s~^ and the spectral index is a = 0.7 (Bicknell et al. 1997). 
Table 1 lists the values of the parameter Vo,3, which provides the scale for the size of the 
velocity function of source size (see eq. M). 

The slopes of both the 6 = 1.5 and 6 = 2 lines are adequate fits to the slope of the 
data points and the overall fits to the data are obviously better for the lower values of tIq. 
The fits could also be improved by increasing the radio power and Thomson optical depth. 
Increasing the radio power increases the brightness temperature of the lobe and hence the 
frequency at which induced Compton scattering becomes important. We have not decreased 
the density still further to produce a more acceptable fit since it is evident from table 1 
that lower densities would involve an unacceptable extrapolation of the fit to the Thomson 
opacity, well beyond the velocities of ~ 1000 km s~^ that the MAPPINGSII code can 
adequately handle. Indeed, the velocity in even the 6 = 1.5, ng = 10 model is greater than 
1000 km s~^ for an overall size <0.5kpc. It can be reasonably assumed, however, that the 
lower densities favoured by our ICS model imply cocoon velocities in excess of 1000 km s~^, 
even if the extrapolation of the shock modeling beyond such a velocity is uncertain. 



4. Discussion 

Our model assumes a uniform distribution of radio emitting plasma in the jet lobes, 
whereas in a number of imaged GPS sources, the emitting region is clumpy; the highest 
brightness temperatures are typically identified with a small number of localized components 
( [Stanghellini et al. 1997| , [Conway et al. 1994| , [Wikinson et al., 1994| ). If these clumps were 



embedded in a lower surface brightness cocoon outside of which the ionizing screen were 
located then the radiation intercepting this screen would be anisotropic. This effect of 
anisotropy is countered however, by the power (typically 1/2.7) entering in the expression 
for Up. Thus, Up would be reduced to approximately 40% of the isotropic value for / ~ 0.1. 
Moreover, as inspection of the images of GPS and CSO sources (the latter all being GPS) 
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reveals, the "clumps" are usually situated at increasing distances from the nucleus with a 
local filling factor of order unity. The beam solid angle corresponding to each clump should 
therefore be of order 2it. Anisotropy need not be an important factor especially in view of 
the insensitivity of the peak frequency to it. 

On the other hand, non-uniformity in the emission increases the brightness temperature. 
In our uniform model, the brightness Temeprature at 5 GHz, is given by: 



\ -2 



'k. (9) 



Six'^xlu'^k V 1027-5 w Hz- V VlOOpc/ V5 x lO^Hz 

and is about an order of magnitude lower than the values observed in subcomponents by 



Stanghellini et al. 19971 , example. A larger brightness temperature of course enhances 



the prospect of induced Compton scattering. 

Interestingly, the Stanghellini et al. study of a complete sample of GPS sources has 
revealed a pronounced peak in the distribution of low-frequency spectral indices near 
a = —1.0, which is in remarkabe agreement with the theoretically predicted spectra 
resulting from induced Compton scattering by external thermal electrons (e.g. |Coppi, 



[Blandford fc Rees 1993| ). This feature, combined with the successful ICS explanation for the 
peak frequency - source size relationship, shows that GPS sources may indeed constitute 
the first real instance where induced Compton scattering plays a relatively clear role in 
extragalactic radio sources. 

We therefore conclude that induced Compton scattering may be an important process 
in the shocked and photoionized plasma surrounding an expanding lobe in GPS sources; 
it can contribute to the formation of the characteristic peak in the radio spectrum when 
the density falls below ~ 10 cm~'^ and free-free absorption becomes inefficient. Expansion 
velocities of 1000 — 2000 km s~^ are then implied when induced scattering is the dominant 
process. Induced scattering may thus be particularly relevant in the high redshift GPS 
sources ( [Elvis et al. 1997| ), which have column depths of the order of a few xlO^^ cm^^^ 



at the lower end of the range estimated by [Bicknell et al. (1997)| on the basis of the FFA 
model. 

In general, the introduction of induced Compton scattering not only increases the range 
of plasma densities over which a low frequency break can be produced, but also increases 
the efficiency of FFA at the lowest frequencies, where photons are being continuously 
replenished as a result of the downscattering of GHz photons. Indeed, the combined effect of 
induced scattering and free-free absorption could also be important in some other sources, 
such as NGC 1275, in which one of the jets appears to be significantly weaker than the 
other at low frequencies. This has been attributed to free-free absorption ( [Levinson, Laor 
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& Vermeulen 1995) and although the observations are on parsec scales in this particular 
source, the brightness temperatures are similar to the values considered here, suggesting 
that induced Compton scattering may well be operating on various jet scales, at the very 
least in conjuction with the more conventional process of free-free absorption. 

We wish to thank Professor Rashid Sunyaev for suggesting that we consider ICS for 
GPS sources and our anonymous referee whose valuable comments and suggestions helped 
to improve the paper significantly. We also thank Chris O'Dea for helpful comments. 
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Figure Caption 



Fig. 1. — A comparison of the predicted and observed relationship between the peak 
frequency and source size for GPS sources. The lines represent the theoretical relationship 
predicted by equation @: sohd and dashed lines are for 6 = 2.0 and S = 1.5, respectively 
(values of no = 1 cm~^ and Uq = 10 cm~^ are indicated) and the other values used are a = 0.7, 
-PsGHz = 10^^'^ W Hz~^ and = lO^^erg s~^ (see the text for definitions of these parameters). 
The data points are from Fanti et al. (1990), O'Dea & Baum (1997)and Stanghellini et al. 
(1997). 
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Size (kpc) 



Table 1: Pcirameters of expanding cocoon 



no S ^0,3 

i 2 1.80 

10 2 0.84 

1 1.5 1.75 

10 1.5 0.81 



